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Abstract 

We study the problem of a random walk on a lattice in which bonds connecting nearest 
neighbor sites open and close randomly in time, a situation often encountered in fluctuating 
media. We present a simple renormalization group technique to solve for the effective diffusive 
behavior at long times. For one- dimensional lattices we obtain better quantitative agreement 
with simulation data than earlier effective medium results. Our technique works in principle in 
any dimension, although the amount of computation required rises with dimensionality of the 
lattice. 
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1. Introduction 

Many naturally occurring diffusive processes do not take place within a static medium. 
Influences of medium fluctuations on transport properties have been encountered particularly 
in polymeric host media (see [6] and references therein). To illustrate such a situation in the 
particular field of molecular biophysics, let us consider the migration of a ligand to a protein 
active site [13], a problem with which we have been concerned recently [9]. Several proteins (e.g. 
myoglobin) have their active site not on the surface, but inside the protein matrix. The ligand, 
in the case of myoglobin a small gas molecule like oxygen or carbon monoxide, has to reach 
its binding site, the heme pocket, by somehow crossing the protein matrix. An analysis of the 
conformational structure of myoglobin [4] suggests that no paths from the outside of the protein 
to the heme pocket are present when a protein is frozen into a static, average conformation. 
However, the protein conformational structure is not static — conformational fluctuations are 
present at physiological temperatures and, in fact, are crucial for the proper functioning of the 
protein [2,4]. There is a good deal of evidence that ligand diffusion cannot take place without 
local volume fluctuations within the protein. If myoglobin within a glycerol- water solvent, for 
example, is cooled to temperatures well below the glass transition of the solvent-protein system 
(~ 200°K), then no ligand diffusion within the protein is observed [1]; presumably, a glass 
transition freezes out conformational fluctuations of the protein, at least on the timescales of 
the experiment, thereby prohibiting diffusion. Furthermore, the ligand almost certainly affects 
the neighborhood through which it is moving. Therefore, a correct treatment of ligand diffusion 
in myoglobin must take into account that local pathways for the ligand will appear and disappear 
randomly. 

We are therefore led to consider the problem of macroscopic transport in a medium in which 
channels are randomly appearing and disappearing over a timescale r. A particular instance 
of such a situation is the case of a random walk that jumps from node to node on a lattice 
where the nearest-neighbor site connectivity fluctuates randomly in time. Simpler still, consider 
a lattice where each nearest-neighbor link opens and closes randomly, or, equivalently, where the 
nearest-neighbor hopping matrix element fluctuates between zero and some nonzero value. This 
last case, although special, is of great interest and wide applicability. A simplified version of 
this process has been used to study ionic conduction in polymeric solid electrolytes and protonic 
diffusion in hydrogen-bonded networks, among other things [6]. It is also clearly relevant to 
ligand diffusion in biomolecules such as hemoglobin [2], although a quantitative theoretical 
study of this situation has not been made. 

For concreteness, first consider the static percolation problem (r = oo) on a lattice in which a 
bond is present with probability p. Macroscopic transport (i.e. the ability of a particle to diffuse 
from one boundary of the system to another) requires that p be greater than pc, the percolation 
threshold of the lattice; = 1 in one dimension, and depends on the lattice structure in higher 
dimensions. When r < oo, macroscopic transport will be possible for any p. An "effective" 
diffusion constant Deff can be defined by Deff = limf^oo(''^(i))/^) provided the limit exists, 
where r(f) is the distance from the origin at time f of a particle that starts at the origin at t = 0; 
Deff will depend on p and r. 
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This problem has been analyzed in previous treatments [6,7,8,12,14], which used effective- 
medium theories for calculating -De//- Qualitatively, their results agree with what one expects 
intuitively: For rapid bond fluctuations (r small), Deff effectively tracks p, while for larger r the 
dynamical problem approaches the static percolation problem, with an increasingly abrupt tran- 
sition from low to high transport coefficient at the appropriate percolation threshold. However, 
although it is well known that effective medium theories can fail quantitatively [5], apparently 
no simulations have been performed up to now to check these results. If effective medium 
approaches fail to give satisfactory quantitative agreement with simulations, a new procedure 
which is capable of doing this is highly desirable, if one ultimately wishes to compare labora- 
tory experiments (say, pressure and temperature studies of ligand diffusion in myoglobin) with 
theory. 

The purpose of this paper is to propose a new treatment of the problem which does precisely 
that. It employs simple renormalization group ideas and is extremely easy to implement on a 
one- dimensional lattice. Furthermore, we compare both this implementation and earlier effective 
medium results with new, detailed numerical simulations. 

2. The Model 

To arrive at a better understanding of the processes discussed in the introduction, we will 
study a simplified problem in discrete space and time. Consider site diffusion on the lattice Z'^. 
At each discrete timestep (which occurs at unit intervals) the particle must attempt to jump 
to a nearest neighbor site. There is a priori an equal probability of l/2c? for the particle to 
attempt to jump to any one of its 2c? neighboring sites. This decision is made independently of 
the state of the bonds. Once the particle attempts to jump in a given direction, the jump will 
be successful if the relevant bond is present ("on"). If the bond is absent ("off"), the particle 
remains at its starting site until the next timestep, when again it attempts to jump in some 
direction. Let hijin) denote the state of the bond from site j to site i at timestep n by hijin) = 
for off and hijin) = 1 for on. Given a particular bond history 



the particle dynamics can easily be cast into the following equation for the time evolution of the 
probability P(i, n \ B) that a random walker is at site i at timestep n: 



Note that the sums are over nearest neighbors. Our case is symmetric in the sense that hijin) = 
hjiin) holds; however, equation (2.2) can be readily used for the case of directed bonds, too. 
This jump process is Markov when the bond configuration is static because all attempts to jump 
are independent of earlier events. However, the jump process is not Markov when the bonds 
fluctuate because the bond state affects the jumps. 



5 = {{&..(0)},---,{M«)},---}, 



(2.1) 




(2.2) 
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We now impose an independent bond dynamics on the bond variables hijin). Each bond 
independently fluctuates on and off at random times as a two-state Markov process that is 
identical for all bonds. The state of each bond at the discrete jumping times is then governed 
by the matrix of one-step transition probabilities Pb'b-, the probability that hijin + 1) = 6' at 
timestep n-\-l given that hijin) = 6 at timestep n. This transition matrix has the general form 

Poo Poi\ _ + {l-p){l-v)\ , . 

where < p < 1 is the probability that any given bond is on and < < 1 is the correla- 
tion factor. Iterating the transition matrix n times then gives the matrix of re-step transition 
probabilities pfjUjin) to be 

Pm{n) Pn{n)\^(pm PnX ^ ( {I - p) + pv"^ {I - p){l - , . 

This shows that when v < 1 temporal correlations decay as z/". The Markov property allows 
us to consider the bond dynamics to be a continuous time Poisson process without loss of 
generality. In this case the correlation time < r < oo (in lattice temporal units) is determined 
by 1/ = exp( — l/r) and the mean lengths of time over which a bond remains either on or off, 
denoted Ton and Tqjj, are given by Ton = ~ p) a^nd Tqjj = t jp respectively. When v = 1 

(r = oo) the transition matrix (2.3) becomes the identity matrix and the bonds become static. 
When 1/ = (r = 0) the bonds are completely de-correlated from one jumping time to the next. 
Because the bond process is Markov, the joint bond/jumping process is also Markov. 

Previous work [7,10] has indicated that at long times the behavior of this and related models 
is diffusive for < p and v < 1 (also see Figure 4 below), namely, it was found that there exists 
a positive constant Dgjj such that the relation 

{r\n))^D,jjn (2.5) 

holds for large values of n. These results strongly suggest, and we will henceforth assume, that 
this effective diffusion coefficient exists and is defined by the limit 



A/f = lim i^L^ = lim 1 V \i\^P{i,n), (2.6) 

where P(i, n) is the probability of finding a particle at the site i at timestep n given that it 
started at the origin at time and that \i\ is the distance of site i from the origin. Thus, P(i, n) 
is the expected value of P{i,n \ B) over all bond histories B, where P{i,n \ B) is computed for 
a given B as the solution of (2.2) that satifies the initial condition P{i,0 \ B) = Here 6{-) 

denotes the Kronecker delta function centered at the origin. It is clear that the P{i,n), and 
hence Dgjj, depend on (p, z^), the parameters associated with the ensemble of bond histories B. 
This dependence will sometimes be indicated by writing P{i,n\p,v) and Dejj{j>,v). While it 
may be that P{i,n) is non-Gaussian [10], we will not address that question here, but rather 
study only the determination of the effective diffusion coefficient Dgjj = Dejj{j>,v). 
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The relation of the Deff defined above in (2.6) to a continuum diffusion description can be 
understood by considering the particle dynamics on macroscopic scales in which the unit lattice 
spacing and unit timestep increment have sizes designated 6x and St respectively. Formally 
taking a macroscopic limit in which 6x and St vanish while {^6xY jSt is held fixed, the macroscopic 
density of particles p = p{t,x) will satisfy the diffusion equation 

The continuum diffusion coefficient is therefore proportional to both Deff and the dimensional 
ratio [6xY/6t, which is invariant under the spatio-temporal scale symmetry (cc ax, t a?t) 
of the equation. 

The value of Deff{j>,v) may be easily determined in the following four limiting cases. First, 
when p = 1 the bonds are always on and the particle dynamics reduces to free diffusion with 
Deff = 1. Because bonds being off can only reduce the transport of particles, it follows that one 
must generally have Deff < 1, with inequality indicating deviations from free diffusion. Second, 
when p = the bonds are always off and clearly Deff = 0. Third, when i/ = the bonds fluctuate 
infinitely fast and the expected state of the bonds is independent from timestep to timestep. 
In that case it is easy to show that Deff = p. Finally, when v = 1 the problem becomes that 
of diffusion on a static, randomly bond-diluted lattice. In this case, the behavior will strongly 
depend on the dimension d of the lattice. For c? = 1 it is easily seen that 

Deff{pA)=\ ^ ' l~ : (2.8) 
otherwise . 

For d > 1 the behavior will be diffusive only \i p > pc, where < 1 is the bond percolation 
threshold for the lattice in question, which depends on d. It is well-known that in this situation 
diffusion is anomalous for intermediate times and the limit (2.6) is approached only for very 
long times; exactly at pc diffusion is anomalous for all times [3,11]. 

In addition to the above limiting cases, it may be readily argued that {j>-,y) ^ Deff{j>,v) 
must be an increasing function of p and a decreasing function of v. However, we would like 
to know quantitatively the full range of behavior of Deff as a function of p and v. We have 
developed a renormalization group scheme that provides a new, non-effective medium method 
for computing this, which will be the subject of the following section. 



3. The Renormalization Group Procedure 

Renormalization group (RG) transformations are best understood as transformations in the 
parameter space of the models in question that are connected with a rescaling of space and/or 
time, and that leave macroscopic properties of the system invariant. For the fluctuating bond 
lattices introduced in the last section the exact macroscopic dynamics would be recovered if 
we could find an RG procedure that coarsens the spatial and temporal units of the lattice so 
that the ratio (^Sxy /St remains fixed while the parameters p and i/ are transformed so that the 
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macroscopic bond fluctuation correlation time rSt and the effective diffusion coefficient Dejj are 
invariant. Of course, we cannot fix Deff exactly because that would require prior knowledge of 
Deff. Rather, what we will do is to match one of the approximates to Deff from the limiting 
relation (2.6) for the fine lattice with the appropriate approximate for the coarse lattice. Upon 
iterating the resulting RG transformation to a fixed point we will obtain an approximation to 
Deff which will prove remarkably accurate in one dimension. 

Let Dn{p,v) denote the n^^ approximate to Deff{p,i') in (2.6), which has the form 

n 



n 



n. ^ — ^ 



(3.1) 



The P{i, n\p,v), and hence Dn{p, i^), can be determined explicitly for any value of n. This task 
is relatively easy when n is small. For example, when n = 1 a, particle can have only either 
moved to a nearest neighbor or remained at the origin. Because the probability of attempting 
to move in a given direction is l/2c? while the probability of the move being successful is p, one 
has 

r (i-p) 
1 



P{i,l\p,iy) 



2d^ 




if Kl = 0, 
if Kl = 1 , 

otherwise , 



(3.2) 



which by (3.1) gives 



(3.3) 



i.i=i 



Similarly, when n = 2 a particle can have only either moved to a next-nearest neighbor, moved 
to a nearest neighbor or remained at the origin. By examining the likelihood of each possible 
path a particle could take to end up at i after two timesteps, one can show 



P{i,2\p,iy) 



j^p^ + {1 - pf + ^p{l - p)u 



1 

2d? 

1 , 



p{l-p){2d-iy) 

,2 



4d^ 




P 



if \i\ 
if \i\ 
if \i\ 
if Kl 



0, 

1, 

2, 



(3.4) 



otherwise , 



which by (3.1) gives 



2d 



p{l -p)v . 



(3.5) 



It is already clear that the complexity of the calculation of the P{i, n\p,v) increases rapidly 
with n. Later we will describe a method for carrying out these calculations for general n. 
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Now we consider a few general properties of Dn- First, each z/) has some of the 

same limiting behaviors enjoyed by Dgjj, namely that Dn is an increasing function oi p and a 
decreasing function of v such that 

Dn{p,Q) = P, i?„(0,z/) = 0, Dn{l,v) = l, Dn{p,u)<p, (3.6) 

with equality only for i/ = or p = 0,1 when n > 1. Second, as i/ is decreased the variation of 
Dn{p,i^) with n becomes smaller and the convergence of (2.6) becomes more rapid. 

We now assume that an evaluation of Dn{p, i^) is given for some fixed n > 1 and introduce an 
TC-step RG transformation of the model parameters (p, i/) to new values (p', i/') that corresponds 
to a coarsening of the spatial and temporal units of the lattice by factors of and n respectively. 
The TC-step RG transformation is chosen to match the bond re-step correlation factor and n^'^ 
approximation to Deff for the original system with the one-step correlation factor and first 
approximation to Deff for the new system. Specifically, this means that by (2.4) 

z/' = z/'^ , (3.7a) 

thereby fixing the macroscopic bond fluctuation correlation time, while by (3.3) 

p' = Di{p',u') = Dn{p,iy). (3.7b) 

It is seen immediately that these equations already represent the RG transformation equations 
for p' and i/' which introduce a flow in the parameter space of (p, z^), see Figure 1. From the 
properties (3.6) of Dn{p, z^), one can easily see that (p, z/) = (0, 1) and (p, z/) = (1, 1) are unstable 
fixed points, and z/ = is a line of stable fixed points under the above RG transformation. 

Our strategy for the approximation of Deff{j>,v) is now straightforward: given re > 1, we 
start with the bare model parameters (p, z/) and iterate the re-step RG transformation (3.7). 
Whenever {jp, v) is not (0, 1) or (1, 1), this process will approach a fixed point of the form {jp* , 0), 
in which case we assign 

D!^f{p,u)^Deff{p\0) = p* . (3.8) 

This value for the effective diffusion coefficient is associated with all models that are connected 
by the same RG trajectory. However, it will be far more accurate for those values of v that 
are near zero, where the convergence of (2.6) is rapid, rather than those near 1. This strategy 
for computing D^'j^{p,v) is carried out in the next section for a one- dimensional lattice and for 
different values of re, and the results are analyzed and compared with simulations. 

4. The RG for a One-Dimensional Lattice 

We now specialize the RG procedure to a one- dimensional lattice, deriving the two-, three-, 
and four-step RG transformations. In the next section we will demonstrate that the almost trivial 
two-step renormalization given below is sufficient to give very good agreement with numerical 
simulations, better in fact than effective medium theories which require considerably more work. 
We will improve upon this agreement by continuing to the three- and four-step renormalization. 
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It should be noted that the number of terms, and hence the labor required to compute the RG 
transformation, increases exponentially with the number of steps used. 

The TC-step renormalization requires knowledge of z/), which is computed from the 

P(^i,n\p,i/) by formula (3.1). For the one- dimensional problem the i runs over E. Fixing {p,i^) 
and noting the general symmetry P(^—i,n) = P{i,n), the probabilities P{i,n) need only be 
computed for i = 1, - ■ ■ ,n. Of course, the probability P(0, n) can be determined from the others 
by the relation 

n 

P(0,7i) = l-2^P(i,7i), (4.1) 

i = l 

but this is not necessary because of its vanishing contribution to (3.1). Below we will use the 
two-step case to illustrate how the P{i, n) can be calculated using a diagrammatic technique to 
classify the various possible paths. 

Although it is not necessary for this calculation, the diagrammatic method introduced here 
for counting particle paths is helpful when more extensive computations are required. The 
diagrams we use are a shorthand for collecting related terms. Consider the diagram shown in 
Figure 2. Each level downward corresponds to a new increment of time. At the top level there 
is only one site shown, indicating that at time m = the particle is at the origin; one level down 
(m = 1) three sites are shown, indicating that the particle can be at sites —1, 0, or 1; and two 
levels down (m = 2) five sites are shown (i = —2 through i = 2, inclusive). Our diagrams then 
trace possible particle paths, with the understanding that at time zero the particle starts at the 
origin. 

All two-step paths are indicated in diagrams (a — /) of Figure 3. In addition to P(2,2) 
and P(l,2), we will compute the probability P(0,2) for illustrative purposes. When the final 
position of the particle is at site i the contribution of each diagram to the calculation of P(i, 2) 
is as follows: 

i = 2: (Diagram a) The particle goes right twice (which has a probabilistic weight of jp^). 

i = 1: (Diagram b ) There are two possibilities for this diagram: the particle first fails to go right 
and then goes right (which has weight — p)pio)', the particle first fails to go left and 
then goes right (weight j[l —p)p). The total weight for this diagram is j[l —p){pio +p)- 

(Diagram c) There are two possibilities for this diagram: the particle first goes right and 
then fails to go left (weight jppoi); the particle first goes right and then fails to go right 
again (weight jp{l — p)). The total weight for this diagram is jp{l — p + Poi)- 

i = 0: (Diagram d) The particle first goes right and then goes left (weight jppn). 

(Diagram e) The particle first goes left and then goes right (weight jppn). 

(Diagram /) There are four possibilities for this diagram: the particle fails to go right 
twice (weight — p)poo)', the particle fails to go left twice (weight — p)poo)', the 
particle first fails to go right and then fails to go left (weight j[l — pY); the particle first 
fails to go left and then fails to go right (weight j[l — pY). The total weight for this 
diagram is i(l - _p)(l - p + poo)- 
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Summing the weights of the appropriate diagrams and eliminating the transition probabilities 
Ph'b using definition (2.3) then gives 

P{l,2\p,u)=\p{l-p){2-u), (4.2) 



which agrees with (3.4) when d = 1. The verification of (4.1) by these probabilities provides a 
useful check of the calculation. Using the probabilities (4.2) to evaluate the second approximation 
to Dejj by (3.1) then yields 

D2{p, u) = P(l, 2\p,u) + 4P(2, 2\p,u) = p- \p{l - p)u . (4.3) 

The approximate effective diffusion constant D^'j^ is then found by iterating the two-step RG 
transformation (3.7). In particular, it is interesting to note that when u = 1 this procedure 
recovers (2.8) exactly. 

A better value of D^f should be obtained when the RG transformation is based on a higher 
order approximation to Dgjj. Applying the above diagrammatic technique to compute all pos- 
sible particle paths, for the three-step case we found 

P(3,3|p,z/)= i/, 

P(2,3|p,z/)= i/(l-p)(3-2z/), (4.4) 
P(l, 3 I z/) = |p(/ + 4(1 - pf) - p{l - p){l - 2p)u - \p\l - py , 

while for the four-step case we found 

P(4,4|p,z/)= 

P(3,4|p,z/)= i/(i_p)(4_3j,)^ 

P(2,4|p,z/)= i/(6- 12p+7/) - |/(2-5p+3/)z/, (4.5) 
P(l,4|p,z/) = ip(4 - 12p+ 15/ - 7/) - |p(4 - 20p + 31/ - 15/)z/ 
- i/ (3 - 7p + Ap^y - \p\l - py - 1/(1 - pfu' . 

The third and fourth approximations to Dgjj are then found by formula (3.1) to be 



D^{p, u) = p- iKi - p)u - \p\i - py , 

Di{p, p) = p- iKi - p)p - \p\i - py - 1/(1 - py - ^/(i - pfp' 



(4.6) 



Because the procedures to calculate these cases are essentially no different from those for the 
two-step case, we omit the details. The effective diffusion constants are again found by iterating 
the three-step and four-step RG transformations (3.7). 
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5. Comparisons with Monte Carlo Simulations 

Extensive Monte Carlo (MC) simulations of random walks in a one- dimensional fluctuating 
bond system were performed. We placed N non-interacting walkers randomly on a chain of L 
sites connected by fluctuating bonds where periodic boundary conditions were employed. At 
each time step the state of each lattice bond was updated in accordance with (2.3), then each 
walker attempted a move to a neighboring site in a random direction. The move was accepted 
if the bond connecting the sites in question was on. This scheme corresponds to the analytic 
model described in Section 2. The simulation was stopped when a walker reached a distance of 
L/2 from its starting point, in order to avoid problems due to the toroidal topology introduced 
by the periodic boundary conditions. For our simulations we usually chose N = L walkers, so 
that each lattice site was occupied with one walker on average. Averaging the square of the 
displacements of these random walkers from their starting positions over M runs provides us 
with averages over walks as well as over bond fluctuations in computing {r^[n)). Typically the 
lattice size L employed was 10"^ and the number of runs M ranged from 10 to 100, depending 
on the quality of the statistics. The results were checked against invariance with respect to 
variations in the parameters L,M, and N . 

A typical result for the mean square displacement [r^in)) is shown in Figure 4. The data 
can be described by 

{T\n)) = D,jjn + Tl{l-cj){n)) (5.1) 

where D^ff is the effective diffusion coefficient, is a residual mean square displacement, and 
4>{n) is a — usually non-exponential — relaxation function with ^(0) = 1 and lim^^oo ^(^) ^ 0. 

In Figure 5 the RG and MC-simulation results for Dgjj are compared; since the small p 
regime cannot be resolved readily when Dgjj is plotted linearly, the data are replotted in Figure 
6 with a logarithmic scale for Dgjj. As can be seen, already the two-step RG transformation 
gives a reasonable agreement with the simulation results, which is improved upon increasing the 
step size. In order to give a more quantitative description of the agreement, we have analyzed 
the relative error {D^j' — D^jp) / D^jp . Figure 7 demonstrates the decrease of the relative error 
with step size for i/ = .999. However, as can be seen, the convergence to zero error is rather 
slow. Surprisingly, the relative error is largest (about 50%) for small values of p, whereas the 
RG gives correct results for p = 0. In addition. Figure 8 demonstrates that the relative error is 
rather independent of i/, particularly for small p. 

It is interesting to compare these results with the prediction of effective medium (EM) 
theories. These theories are usually expressed in terms of the correlation time r, which is 
related to i/ by 

i = log(i). (5.2) 

We note that various EM approaches give somewhat different predictions; e.g. [14] predicts a 
scaling Dgff oc p^ /t for small p and large r, whereas [7] predicts Dgff oc p/r, in agreement with 
our MC-simulations and RG results. Therefore we choose to compare the results of [7] with our 
simulations. The one- dimensional self-consistency equation of [7] can be evaluated analytically 
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with the result 

= 1 + 2r(l - pf - ^{l + 2T{l-pff -p{2-p). (5.3) 

Note that this equation gives the correct results for p = 0, 1 and for r = (alternatively, for 
z/ = by (5.2)). Figure 9 shows the relative error of this prediction for Dgjj with respect to 
our simulations. Notice that EM theory systematically underestimates the effective diffusion 
coefficient, the error increasing with r, particularly in the intermediate to large p regime, i.e., 
below the percolation threshold in one dimension. Everywhere, the relative error is larger even 
than that of the two-step RG results. 

6. Discussion 

We have presented a renormalization group approach to obtaining an approximate effective 
diffusion coefficient for random walks on a fluctuating lattice. This procedure was applied to a 
one- dimensional lattice, where it is relatively simple to implement, and was found to be in good 
quantitative agreement with Monte Carlo simulations. The results can be improved by taking 
into account a larger step size in the renormalization procedure. It might be hoped that, e.g., 
using a multi-step renormalization scheme, the procedure presented could be modified in order 
to account also for nondiffusive effects; see (5.1). 

The application of this approach to higher dimensions — two is already interesting — is 
straightforward, but requires much more calculation than the one- dimensional problem consid- 
ered above. In addition, there appears to arise a problem of principle: It is reasonable to expect 
that, as 1/ ^ 1, one should be able to recover the percolation limit. Specifically, as z/ ^ 1 while 
holding p fixed, we expect the limiting diffusion coefficient as a function of p should vanish for 
p < Pc, and be positive for p > pc- This behavior should be reflected in a renormalization group 
flow like it is sketched in Figure 10. In particular, there should arise a nontrivial fixed point at 
V = 1 and p = Pc- However, from the general properties (3.6) of the Dn one can immediately 
conclude that such a fixed point cannot arise for any finite n in our renormalization procedure. 
So, the best that can be hoped for is that the proper behavior is approached as the number of 
steps in the renormalization goes to infinity. Said another way — because of the dimensional 
dependence of pc, one must take enough steps in the renormalization procedure to "see" the 
dimensionality of the lattice. The two-step procedure, employed so successfully in one dimen- 
sion, does not give a very good approximation in two dimensions. However, the analogue of the 
four-step procedure, which is the minimum needed to see simple closed paths, begins to see the 
proper trend. 
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Figure Captions 

1: Renormalization group flow for d = 1; crosses denote fixed points. 
2: A nodal tree for n = 2. 

3: Diagams a — f showing all the possible paths ending at i = 2,1,0. 

4: Average mean square displacement {r^{n)) vs n; parameters as indicated; the dotted line 

denotes the long-time behavior, see (5.1). 
5: Comparison of 2-step (alternating lines), 3-step (dashed lines), and 4-step (solid lines) RG 

results, and MC-simulations for Deff vs p for i/ = 0, 0.9 (circles), 0.99 (squares), 0.999 

(triangles). 

6: Same as Fig. 5, with Deff on a logarithmic scale in order to resolve the small p regime. 
7: Comparison of the relative error of the 2, 3, and 4-step RG results vs p; i/ = 0.999. 
8: Relative error vs p of the 4-step RG results for u = 0.9, 0.99, 0.999. 
9: Relative error vs p of effective medium results for i/ = 0.9, 0.99, 0.999. 
10: Renormalization group flow for d > 1; crosses denote fixed points. 



Figure 2: A nodal tree for 
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Figure 3: Diagrams a — f showing all the possible paths ending at i = 2,1,0. 



